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1.0  INTRODUCTION 

The  program  designs  manufacturing-quality  control  systems  to  mini- 
mize the  expected  operating  cost  of  a single  product  for  a user  speci- 
fied Average  Outgoing  Quality  Limit,  (AOQL).  The  costs  considered  are: 
raw  materials,  inspection,  repair  or  replenishment  of  defectives,  man- 
ufacturing and  assembly.  The  sampling  plans  are  single  sampling  by 
attributes  on  a lot-by-lot  basis.  The  lot  size  is  specified  by  the  user. 

The  program  is  written  in  FORTRAN  for  the  UN  I VAC  1108  at  the  Illinois 
Institute  of  Technology  and  was  designed  for  ease  of  modification  for 
other  computer  installations.  The  program  uses  forward  recursion  dynamic 
programming  to  solve  the  problem.  The  decision  variables  include  the 
sampling  plans  (n  and  c),  the  fraction  defective  generated  oy  a manu- 
facturing or  assembly  process,  and  the  quality  characteristics  of  the 
raw  materials. 

The  underlying  cost  models  assume  that  a large  number  of  lots  of 
the  products  are  to  be  processed  and  the  decision  rules  are  independent 
of  the  actual  outcome  of  the  rules  at  previous  stages.  The  generation 
of  additional  defects  at  a stage  is  assumed  to  be  independent  of  the 
defects  generated  at  previous  stages.  Inspection  is  assumed  to  be  per- 
fect, and  defects  generated  at  any  previous  stages  are  also  detected 
if  a unit  is  inspected. 

Section  2 describes  the  problem,  solution  approach  and  terminology 
used.  Section  3 is  an  overview  of  the  computer  program.  Section  4 
specifies  the  input  format.  Section  5 describes  the  output  and  inter- 
pretation. Section  6 describes  possible  modification  of  the  program. 
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2.0  PROBLEM  DESCRIPTION  AND  TERMINOLOGY 


! 

The  model  may  be  composed  of  the  following  types  of  stations: 

Station  Type 
Inspection 
Manufacturing 
Assembly 

In  this  section  each  of  the  station  type's  cTJ^t  will  be  defined 
and  the  recursive  relationship  shown.  First  the  basic  solution  approach 
will  be  discussed  in  2.1 . 

2.1  The  General  Recursive  Relationships 

The  program  uses  forward  recursion  dynamic  programming.  The  state 
variable  linking  the  successive  stations  is  the  average  fraction  defec- 
tive in  the  lot.  To  discuss  a general  station,  let  us  call  the  average 
incoming  quality  the  pj  and  the  average  outgoing  quality  AOQ.  The 
minimum  cost  for  all  the  stations  up  to  and  including  the  one  under  con- 
sideration, say  n,  is  denoted  by  fR (AOQ ) and  depends  upon  the  value  of 
AOQ  resulting.  The  basic  recursive  relation  is 

fn(A0Q)  = min  / stage  n cost  + f ,(  p,)  \ 

n decision  * n-l  i J 

variables 
of  station  n 

where  AOQ  is  a function  of  pj  and  the  decision  variable  at  that  station. 
The  program  implementation  establishes  a grid  of  values  on  the  state 
variables  and  any  continuous  decision  variables. 

2.2  Inspection  Station 

An  inspection  station  sampling  plan  is  specified  by  (n,c).  A random 
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sampling  of  size  n is  selected  from  the  lot  (of  size  N).  If  the  number 
of  defectives  in  the  sample  is  c or  less,  the  discovered  defects  are 
repaired  (or  replaced  with  good  ones)  and  no  more  inspection  is  done 
on  that  lot.  If  there  are  more  than  c defectives,  the  entire  unsampled 
portion  of  the  lot  (N-n)  is  inspected  and  all  defectives  repaired.  This 
is  called  a sinqle  sampling  plan  with  rectification;  the  quality  of  per- 
formance will  be  at  least  as  good  after  sampling  as  before.  Symbolic- 
ally, the  inspection  station  may  be  represented  as  in  Figure  2.1. 


PI  I 

n,c 

f(Pj) 


AOQ 

f(AOQ) 


Figure  2.1  Inspection  Station 


For  a particular  value  of  AOQ,  the  utility  of  the  inspection  is 
considered  by  comparison  of  fn(A0Q)  with  fn_^(A0Q).  Note  that  for  an 
inspection  station  to  have  a feasible  plan,  AOQ  < Pj 

For  the  final  inspection  station,  the  average  outgoing  quality 
limit  AOQL  is  used  instead  of  AOQ  as  the  state  variable  and  fn(A0QL)  is 
the  minimum  average  total  cost  associated  with  the  AOQL.  The  plan  (n,c) 
for  the  final  inspection  station  will  guarantee  that  the  average  outgoing 
quality  AOQ  for  any  value  of  Pj  is  less  than  or  equal  to  AOQL. 

2.2.1  Inspection  Station  Costs 

The  expected  cost  of  inspection  at  an  inspection  station  is  denoted 
by  kg  • ATI  where  ATI  is  the  average  total  inspection:  ATI  = n + 
(N-n)(l-P  ) where  P is  the  probability  of  accepting  the  lot.  The  cost 

d a 
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per  unit  inspected  is  denoted  by  k^.  The  expected  repair  cost  is 
denoted  by  k^  • ATI  • Pj  where  Pj  is  the  fraction  defective  in  the  lot 
and  hence  ATI  • Pj  is  the  expected  defectives  to  be  repaired,  k^  is  the 
cost  per  unit  repaired.  The  expected  cost  at  an  inspection  station  is 
of  the  form  (k^  • Pj  + k^jATI.  The  probability  of  acceptance  is 

z ("Pj)  e_npI 

x! 

x--0 

2.2.2  Inspection  Station  Recursions 

The  inspection  station  could  be  associated  with  the  inspection  of 
raw  materials  or  after  a manufacturing  or  assembly  station.  If  the 
inspection  is  of  a raw  material  and  no  manufacturing  has  been  done  on 
the  component,  then  the  recursive  relation  is 

f^AOQ)  = min  {k^N  + (k^  + k2)  ATI} 

Prn,c 

where  k^  is  the  cost/unit  of  the  raw  materials  and  may  be  a monotoni- 
cally  nonincreasing  function  of  pj.  The  resulting  Average  Outgoing 
Quality,  AOQ,  equals  P,P,(N-n/N). 

If  the  inspection  is  of  a component  which  has  one  or  more  previous 
manufacturing  or  assembly  stages,  the  recursive  relationship  is 

f n ( AOQ ) = min  {(k^  + k?)  ATI  + f^^Pj)} 
n,c 

2.3  Manufacturing  Station 

At  a manufacturing  station,  each  of  the  N units  in  the  lot  are 
manufactured.  Symbolically,  the  manufacturing  station  may  be  represented 


4 


as  in  Figure  2.2.  The  manufacturing  operation  generates  additional 
defectives  with  a probability  of  p1  for  each  unit.  Thus,  the  resulting 
output  has  a quality  of  conformance  n<_  better  than  that  of  the  input, 
and  probably  worse. 


f(pr)  H PT) 


Figure  2.2  Manufacturing  Station 


2.3.1  Manufacturing  Station  Cost 

If  the  cost/unit  manufactured  is  the  manufacturing  cost  for  the 
whole  lot  is  k^N.  The  cost/unit  k3  may  vary  as  p'  is  changed.  Assume 
that  the  value  of  k^  is  monotonically  non-increasing  with  p'  (i.e.  k^ 
stays  the  same  or  decreases  as  p'  increases). 

2.3.2  Manufacturing  Station  Recursions 

Let  fn  -j  (p j ) be  the  minimum  expected  cost  for  the  station  previous 
to  stage  n as  a function  of  the  fraction  defective  incoming  to  stage 
n,  Pj.  Then  the  recursive  relationship  is 

fn(PT)  = min  {N  ’ k3  + fn-l^pI^} 

p' 

where  p-j.  is  the  average  outgoing  quality  level  for  a manufacturing  stage 

PT  = p*  + Pj  - p'pr 
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2.4  Assembly  Station 


At  an  assembly  station  two  or  more  components  are  assembled  (and/or 
processed)  into  a single  component.  The  aggregation  of  the  components 
results  in  an  Aggregate  Incoming  Quality,  pAJ.  Furthermore,  the  assembly 
and/or  processing  of  the  components  may  generate  additional  defectives 
with  a probability  of  p'.  Figure  2.3  represents  this  symbolically. 


Figure  2.3  Assembly  Station 

The  quality  of  incoming  component  x is  called  pT  . 

l ,x 

2.4.1  Assembly  Station  Costs 

Let  by  the  cost/unit  of  the  assembly  operation.  Then  the  assem- 
bly cost  is  Nkg  because  the  entire  lot  must  be  assembled.  This  cost/ 
unit  may  be  allowed  to  vary  as  p'  (the  probability  of  a defective  gen- 
erated by  the  assembly)  varies.  Assume  that  the  value  of  k^  is  monotoni- 
cally  non-increasing  function  of  p'. 

2.4.2  Assembly  Station  Recursions 

The  assembly  station  recursion  consists  of  two  phases  for  each 
assembly  station.  In  the  first  phase,  f(pAI)  Is  determined.  In  the 
second  phase,  the  recursion  is 


fn(PT)  = min  lNkb+f(pAI)} 

where 

PT  = P*  + PAI  - P'PAI. 

In  the  first  phase,  the  components  are  aggregated  first  by  aggre- 
gating the  first  two  components  to  obtain  an  aggregated  f(pAj).  Then 
the  result  is  aggregated  with  the  third  component  to  get  a new  T(pAj)- 
This  continues  aggregating  one  additional  each  time  until  all  are  aggre- 
gated, and  the  result  is  then  called  f(pAj).  The  general  step  in  this 
phase  would  be 


f(pAI)  = min  <f(Pi,x)  + T(pai)> 


. f0RTPAK  V on  the  «W«  «»  * "11n0iS 
The  prog™  -*  Ja5  desi9„ed  for  ease  of  use  at  f-M- 

4 <«  Section  4.1.  the  multi -component, 

cribed  in  Sectioi  program  for  tne  mu 

The  manual  describes  the  us  ^ ^ i5  one  or  more  assembly 

multi-stage  manufacturing  syste  ^ Mmlfacturing  system  is  thus 

stages-  The  singie-coWoneo« ■ » ^ lt  is  assumed  that 

each  station  has  at  mo 

have  mre  thah  one  direct  ^ (perhap5,  partial  - 

The  term  coaonent  is  used  ^ ^ the  output  (component) 

pleted  product  resulting  from  any  s ^ ^ direct  successor. 

for  one  station  is  the  input  ^ #f  dy„amic  programing-  it 

The  tem  stage  is  from  t e chart  (section  4.2). 

refers  to  the  column  location  on  a ^ ^ „Mch  is  descrihed 

The  program  has  an  elahora  structure  are  surname 

in  section  4.0.  The  basic  rules  ^ ^ ? ^ ^ programing 
in  3.1.  The  algorithms  descr,  e ^ the  algorithm  for  a manu  a 

forward  recursions.  called  to  find  f(PM> 

- rvr  • - - - ~ - 

assembly  component 

determine  the  minimum  cost  for  PT^  ^ stations  except  for  the 

Since  the  ATI  term  is  use  a ,.t  of  tables  which  resu 

fi„al  inspection,  it  is  convenient^  ^ ^ ^ Itisp,  does 

1n  minimum  ATI  for  various  va  ues  ^ AQQ  are  done  by  sub 
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Upon  processing  of  the  station  for  the  final  product,  backtracking 
must  be  done  through  the  sequentially  generated  station  solution  table 
in  order  to  find  the  optimal  solution.  This  is  done  by  the  subroutine 
BAKTRK. 

3.1  Consistency  of  Assembly  Chart 

The  following  conditions  must  be  met  when  using  the  program.  They 
are  referred  to  as  the  rules  for  the  consistency  of  the  assembly  chart. 

1.  There  must  be  at  least  two  stages. 

2.  The  numbering  of  stages  and  components  must  be  consecutive 
starting  from  one,  i.e.  1,  2,  3,...,  etc. 

3.  There  must  be  at  least  one  station  for  each  component  number 
and  for  each  stage  number. 

4.  No  two  Inspection  stations  can  be  placed  consecutively  along 
a component. 

5.  For  an  assembly  station  at  least  two  components  must  be  assem- 
bled. 

6.  There  can  be  one  and  only  one  component  resulting  from  an 
assembly  station  and  for  the  resulting  component  there  must  be  at  least 
one  station  (of  any  type)  in  one  of  the  following  stages.  The  resulting 
component  number  must  be  one  of  the  incoming  component  numbers  for  the 
assembly  station  under  consideration. 

7.  The  last  stage  of  the  manufacturing  system  must  be  an  inspec- 
tion station. 

3.2  State  and  Decision  Variable  Grids  and  FMAX  Value 

Since  the  algorithm  employs  dynamic  programming  and  the  state  vari- 
ables and  decision  variables  (other  than  n and  c)  are  continuous,  it  is 
necessary  to  provide  a finite  grid  of  values  for  the  continuous  variable. 
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In  practice  it  makes  little  sense  to  use  a grid  smaller  than  0.5% 
because  of  the  accuracy  of  controls  and  cost  estimates.  The  user  must 
provide  the  grid  size  desired  for  these  variables.  The  same  grid  size 
is  used  for  all  state  variables.  The  grid  sizes  for  the  decision  vari- 
able for  manufacturing  and  assembly  stations  may  be  different.  The  pro- 
gram uses  0.5%  for  the  finest  possible  grid  size  and  rounds  all  values 
on  that  basis. 

The  maximum  limit  on  the  average  fraction  defective  from  a lot  is 
conceptually  100%.  Usually  the  range  of  interest  is  much  smaller.  The 
program  allows  the  user  to  specify  a maximum  which  is  smaller  than  100%. 
This  allows  a possible  suboptimal  solution  to  be  obtained.  The  limit 
is  specified  by  employing  a constant  multiplicative  factor,  called  FMAX 
on  the  value  of  AOQL.  For  example,  if  the  AOQL  value  is  0.05  and  FMAX 
value  is  2.0,  then,  the  maximum  value  of  the  state  variable  is  0.10 
(FMAX  * AOQL)  throughout  the  system.  It  should  be  noted  that  the  larger 
the  FMAX  value  is,  the  more  unlikely  a suboptimum  becomes.  However 
the  price  is  that  the  computation  time  is  prolonged. 

3.3  Explanation  of  the  Program  Elements 

The  program  has  23  elements,  one  of  which  is  the  MAIN  program. 

Every  statement  in  each  element  has  a sequence  of  numbers  for  ease  of 
ordering.  The  first  three  characters  of  the  statement  number  represent 
an  abbreviation  of  the  name  of  the  element.  The  following  table  gives 
the  names  of  the  elements  (subroutines),  where  they  are  called  from,  their 
abbreviations,  number  of  statements  in  each  element  and  their  basic  func- 
tions. 
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Table  3.1  Program  Elements 


Program 

Element 


; MAIN 


C0NTRL 


Called 

From 


Abbre- 

viation 


# of 

Statements 


INPUT 


ASSEMB  INPUT 


C0STFN  PR0CES 

ASSEMB 
RAWMAT 

RANK  PR0CES 

RAWMAT 
ASSEMB 


Algorithm  Subroutines 


Basic  Function  Performed 


Sets  parameters  and  calls 
other  subroutines. 


Controls  calling  of  sub- 
routines according  to  the 
structure  of  a problem. 

It  also  defines  the  Ran- 
dom Access  File  and  con- 
trols output. 


Reads  input  and  controls 
other  input  reading  sub- 
routines. 


Reads  inputs  costs  for 
raw  materials. 


Reads  input  costs  for 
inspection  station. 


Reads  input  cost  for 
manufacturing  station. 

Reads  input  costs  for 
assembly  station. 

Generates  cost  coordinates 
according  to  the  specified 
cost  functions  1 ,2 9. 

Ranks  the  set  of  cost  co- 
ordinate pairs  to  deter- 
mine whether  or  not  they 
are  monotonically  non- 
increasing function  of  p. 


P0ISSN 

INSP1 

PSN 

17 

Calculates  the  probabil- 

ity  of  acceptance  based 

on  Poisson  Distribution. 

INSP1 

AT  INC 

IP1 

73 

Generates  tables  of  ATI, 
n and  c. 

I 


1 


INSP2 

C0NTRL 

IP2 

125 

Algorithm  for  inspection 
station. 

INSFIN 

C0NTRL 

IFN 

71 

Algorithm  for  Final  In- 
spection Station. 

ASSMBY 

C0NTRL 

ABY 

68 

Algorithm  for  assembly 
station. 

PR0DN 

C0NTRL 

PDN 

87 

Algorithm  for  manufactur- 
ing and  assembly  station. 

Subroutines  for  Output 

DISK 

C0NTRL 

DSK 

39 

Writes  station  outputs  to 
the  Random  Access  File. 

BAKTRK 

C0NTRL 

BTK 

93 

Backtracking  according  to 
the  structure  of  the  pro- 
blem obtains  and  reports 
optimal  solution  from  the 

Random  Access  File. 

DATOUT 

MAIN 

DAT 

75 

Prints  some  outputs. 

Auxiliary  Subroutines 


HEAD 

C0NTRL 

DAT0UT 

INPUT 

HED 

10 

Prints  heading. 

SETMX 

C0NTRL 

SMX 

19 

Arranges  data  before 
calling  INSP2. 

AT  INC 

MAIN 

ANC 

82 

Arranges  data  and  calls 
INSP1 . 

ZER0IN 

INSP1 

ZEN 

52 

Searches  for  value  of  n 
for  each  values  of  AIQ, 
AOQ,  and  c. 

FCT 

ZER0IN 

FCT 

72 

Function  subroutine. 

3.4  Size  Limitations  of  Program 

There  is  really  no  limit  to  the  size  of  the  problem  that  the  program 
can  solve  in  terms  of  the  number  of  components  and/or  stages  unless  the 
capacity  of  the  computer  is  exceeded. 
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However,  it  should  be  pointed  out  that  a maximum  of  only  999  com- 
ponents of  each  type  can  be  accommodated  in  any  stage.  Since  there  is 
no  limit  on  the  number  of  stages,  we  can  always  create  additional  stage(s) 
to  accomnodate  any  number  of  excess  stations  of  one  type  beyond  999. 

For  example,  if  we  have  a problem  which  requires  1500  components 
and  each  of  these  components  has  to  be  inspected  as  soon  as  it  enters 
our  system.  The  problem  may  look  like  Figure  3.1. 


Stage 

Component 

1 

2 


1500 


O 

□4- 


O 


Figure  3.1  Wrong  Arrangement  for  Size  Violation 


In  order  to  satisfy  the  constraint,  we  can  create  an  additional 
stage  and  perhaps  arrange  the  stations  as  in  Figure  3.2. 


Stage 

Component 

1 

2 


999 

1000 

1001 


o- 

O- 


of 


Of- 

O 


1500 


Of- 


Figure  3.2  Correct  Arrangement  for  Size  Violation 


13 


4.0  INPUT 


» 


i 


This  program  package  consists  of  23  elements,  2 of  which  can  be 
modified  by  the  user  in  order  to  accomodate  a bigger  size  problem  (e.g. 
more  components  or  stages)  or  to  fit  their  own  type  of  cost  function. 

The  two  elements  are 

1.  MAIN  Program 

2.  SUBROUTINE  C0STFN  . 

See  Section  6:  Possible  Modifications  of  the  Problem. 

The  users  are  required  to  prepare  input  data  and  place  them  on  a 
sequential  data  file.  The  program  reads  the  input  data  and  provides 
diagnostic  messages  (if  any)  so  that  the  user  can  correct  their  input 
easily.  Diagnostics  are  usually  the  result  of 

(i)  miss-typing  of  input  statements. 

(ii)  dimensions  of  some  arrays  exceeded. 

If  there  is  any  error  that  the  program  can  detect,  the  program 
will  not  proceed  to  calculate  the  optimal  solution  but  instead  a partial 
output  with  diagnostic  messages  will  be  provided  to  identify  the  sources 
of  the  errors. 

4.1  Schematic  Diagram  of  Input/Output  Files 
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IMIi 
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The  Input  Data  File  is  assumed  to  be  sequential  and  is  read  by  the 
1/0-variable-name  called  NRD.  Each  record  is  assumed  to  be  of  80  bytes 
long  (similar  to  an  80  column  card).  The  file  is  formatted. 

The  Random  Access  File  is  created  by  the  Fortran  DEFINE  FILE  state- 
ment. The  program  will  determine  the  size  of  this  file.  Its  1/0  vari- 
able is  NDS.  If  the  program  detects  any  error  (hence  computation  will 
be  terminated  by  the  program)  this  file  will  not  be  created.  It  is  to 
be  created  as  an  unformatted  48-byte  record  file,  if  there  is  no  error. 

The  Output  File  assumes  the  standard  133  character  record  file  for 
the  paper  printouts.  Its  1/0  variable  name  is  NWR. 

4.2  To  Create  Input  for  a Problem 

Since  there  are  only  3 types  of  stations  allowed  in  a problem 

i.e.  Manufacturing,  Assembly  and  Inspection.  The  letters  M,  A,  I, 
represent  them,  respectively. 

Every  station  must  be  identified  by  component-number  and  stage- 
number  as  well  as  its  type.  In  addition  a station  may  or  may  not  have 
a name  associated  with  it  depending  on  its  type.  The  following  rules 
must  be  adhered  to: 

1.  Each  station  in  the  system  is  allowed  to  have  a name  with  a 
maximum  of  4 characters. 

2.  Every  assembly  station  must  have  a name.  This  is  necessary 
since  two  or  more  components  will  be  assembled  into  one  and  hence  it  is 
necessary  to  know  which  station  assembles  them. 

3.  Naming  for  M and  I stations  is  optional  because  such  stations 
are  uniquely  identified  by  their  component- numbers  and  stage-numbers. 

Consider  a simple  manufacutring  system  in  which  two  components 
are  assembled  into  one  unit  as  in  Figure  4.2. 
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Component  1 


Component  2 

Stage 


AL 


1 


DAN 


Figure  4.2  Example  Assembly  Chart 


The  assembly  chart  is  constrained  to  be  as  compact  as  possible  to 
the  left  and  the  components  and  stages  are  numbered  accordingly.  This 
manufacturing  system  can  be  described  as  having  two  components  and  four 
stages.  In  the  third  stage,  the  first  component  after  manufacturing 
is  assembled  with  component  2,  and  the  assembled  unit  is  inspected  at 
the  inspection  station  called  MIKE.  This  component  could  be  numbered 
either  1 or  2 - See  Section  3.1.  (See  also  STRUCTURE  Section  4.3.2.) 

We  identify  the  problem  as  in  Table  4.1. 


Table  4.1  Example  Structure  Naming 


Component 

Stage 

Station  Type 

Station  Name 

1 

1 

I 

ALW»* 

1 

2 

M 

DAN* 

1 

3 

A 

JM 

2 

1 

I 

»m 

2 

3 

A 

3m 

2 

4 

I 

MIKE 

* V = blank 


Suppose  that,  in  reality,  we  are  producting  a lot  size  of  1000  units 
and  two  units  of  component  #1  are  required  for  each  assembly.  Since  all 
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lot  sizes  must  be  the  same,  then  we  create  another  ficticious  component 
called  #3  which  is  the  same  item  as  #1,  each  with  a lot  size  of  1000. 
Figure  4.3  shows  how  the  assembly  chart  would  have  to  be  modified. 


Component  1 
Component  2 
Component  3 
Stage 


AL  DAN 

12  3 4 


Figure  4.3  Example  Assembly  Chart  - Modified 


Although  in  Stage  1,  we  have  2 inspection  stations  called  "ALOO" 
they  are  not  taken  as  the  same  station  by  the  program  because  they  have 
different  component-numbers , i.e.  1 and  3.  Similarly  for  the  process- 
ing stations  in  stage  2,  they  are  treated  independently  even  though  they 
have  the  same  name,  DAN.  Therefore  cost  data  must  be  given  for  each  of 
them  seperately,  but  the  data  will  be  the  same. 


4.3  Input  Sections 

The  input  data  required  by  the  program  may  be  divided  into  sections 
as  follows: 

1.  QUALITY  CHARACTERISTICS 

2.  STRUCTURE 

3.  INSPECTION  COST 

4.  RAW  MATERIAL  COST 

5.  MANUFACTURING  C0ST 

6.  ASSEMBLY  C0ST 

It  is  instructive  to  consider  each  record  in  the  sequential  input 
data  file  as  an  80-column  card. 
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Note:  Columns  61-80  are  reserved  for  comments. 


Each  input  section  has  a title  card  and  the  last  card  in  the  data 
file  is  the  ENDATA  card  as  shown  below. 


Each  section  will  now  be  discussed  in  detail. 

4.3.1  Quality  Characteristics 

This  section  must  be  first  in  the  data  file.  It  consists  of  two 
cards  only.  The  first  card  is  the  title  card.  Ks  content  is: 


Zalum 

1 2 3 4 S 6 7 8 

l234567R90123'6G7r.901234  567  C901234  567C901?3':66709;>12::45C78<!012345078''ri:'3o07i;V(1 

QUALITY  CHARACTERISTICS 

ComiKiits 

The  second  card  format  is: 


Qolunn 

1 2 3 4 5 6 7 8 

1234 567890 12 34 5678901 2 34 5678901 234 5670901 2 345678901 2 34 5678901 2 34 5678901 234567890 

SX 

F5.1 

5X 

FS.3 

SX 

TS.3 

lx 

18 

« 

b 

— 

c 

d 

Comments 

a = Final  A0QL  desired  for  the  problem 
b = Grid  size  to  be  used  for  state  variable 

c = FMAX,  this  is  a factor  which  establishes  the  range  of  p values  which 
the  program  will  examine  to  find  the  optimal  solution.  (See  Section 
3.2) 

d = The  lot  size  of  the  final  product 
4.3.2  Structure 


This  section  describes  the  manufacturing  stations  involved  in  the 
system  to  the  computer.  Although  it  is  natural  to  input  stations  in  the 
order  of  components  or  stages,  the  program  can  handle  any  order  the  user 
may  desire. 

The  general  input  format  for  the  STRUCTURE  input  section  is: 

Title  card  


;dur  i 

1 2 3 4 5 6 7 C 

l?346r7.*o,- >l2M5S7f."0!  237S.'7r.(>(il23--  5578901  :6f.70fO;?:-V.f.'r.'!'Ji:345S'»M'  u.TL'‘f. 

STRUCTURE 

Com  n;  r.ts 
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I I ^ 1 Whti 


Subsequent  Cards 


> 


Column 

1 2 3 4 5 6 7 8 

I2345C  70o01234r'C7890173‘'6*'7Co01734rj6789012345678*O7?34S(.~f.*101?345G73o{'1734  -.678*0 

SI 

15 

15 

3 

wm 

b 

fl; 

ll 

A4 

* 

fl 

B 

B 

1 

b: 

» 

a: 

i 

b: 

D 

ID 

Coflmmts 

where  a = blanks 

b = component  number 
c = stage  number 
d = type  of  station  (M.A.I) 
e = station  name  (Alpha  Numeric) 


Each  card  specifies  a component,  therefore  it  can  contain  a maximum 
of  5 stations  belonging  to  the  same  component.  For  the  example  of  section 
4.2,  the  complete  input  data  is  as  shown  below: 


Column 

1 2 3 4 5 6 7 8 

123456789012345678901 2345678901 2345678901 23456789012345678901 2345678901234567890 

STRUCTURE 

1 HALM  2HOAN6  3AJ0ER  FIRST  COMPONENT  STNS 

2 1IHW  34 JOE R 4 MICE  SECOND  COMP  STATIONS 


If  there  are  more  than  5 stations  for  a component,  repeat  the  format 
until  all  the  stations  along  the  component  are  completed.  If  you  omit 
the  component  number,  the  program  wi 11  assume  that  the  station  belongs  to 
the  last  component  it  has  read. 

4.3.3  Inspection  Costs 

These  are  units  costs  associated  with  an  inspection  Station.  Costs 
associated  with  all  inspection  stations  must  be  grouped  in  this  section. 
Each  card  represents  an  inspection  station. 


i 

i 


r 

\ 

h 

f 

i 
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Title  Card 


bl«i 


I 


- » 


itsrccTiON  costs 

funit: 

Subsequent  Card 


tolu 


M 

IS 

IS 

SI 

no.  2 

no. i 

• 

ft 

c 

4 

r—  — 

where 


a * component  number 
b = stage  number 
c » unit  repair  cost 
d = unit  inspection  cost 


4.3.4  Raw  Material  Costs 

The  basic  components  entering  the  manufacturing  system  are  called 
raw  materials.  Their  unit  costs  corresponding  to  qualities  (p)  are 
entered  in  this  section.  See  descriptions  of  the  cost  functions  in 
Section  4.4. 

The  general  format  is 
Title  Card 


i ' t'  "l  ft  l • » « 

l7MS67MQ1?34GG7n90173456789012  JftS470S012Jft  SftTOTOlJMSC 78901 2Sft567ftft01?M547W0 

M MTtllAL  COSTS 

Subsequent  Cards 

Each  component  of  the  raw  materials  requires  at  least  2 cards  to 
specify  its  coordinates.  The  first  card  identifies  component  number, 
types  of  cost  functions  and  number  of  cost  coordinates.  Its  format  Is 
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i 


fcolmm 


12  3 4 5 6 7 8 

1234567f f ft ) 7 Kl, (.  7r<1Q1, •'34!,ri7890i:>34  567390)7.17  78"0 1734  '■C/WWIMlVU 


sx 

IS 

IS 

IS 

• 

b 

C 

Cooncnts 

where  a * component  number 

b * type  of  cost  function  (0,1,..., 9) 
c * number  of  coordinates  if  the  cost  function  Is  of  type 
0;  otherwise  leave  blank. 

Other  cards  following  this  card  have  the  format: 


fcoluin 


1 2 3 4 5 6 7 C 

h ? 34 r 7C°P1  ? 3*r.C.7r«01  ? V rjf  7B»;il 7 y 1 . ..'rr.7'inpr3,fM  V.<Kli:3'i5('TP.r  0173' 


iox 

F5.3 

FS.J 

FS.J 

FS.J 

FS.J 

F5.J 

F5.J 

F5.  3 

F5.3 

FS.J 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

Contriv'd  ts 

*See  Section  4.4  for  explanation  of  these  fields. 

4.3.5  Manufacturing  Costs 

Cost  associated  with  the  manufacturing  stations  are  Inputted 
through  this  section.  The  form  of  Input  Is  similar  to  that  of  RAW 
MATERIAL  C0STS  except  that,  in  this  case,  we  have  to  identify  the  com- 
ponent and  the  stage  of  the  station. 

Title  Card 


lolumn 

1 2 3 4 S 6 7 8 

1234567890123456789012345678901234567890123456789012345678901 2345678901 234567890 

IWNUFAC TURING  COSTS 

Coanenti 

where 


a = component  number 
b = stage  number 
c = type  of  cost  function 

d = number  of  cost  coordinates  If  c = 0;  otherwise  leave  blank 
Other  cards  have  the  following  format: 


♦See  Section  4.4  for  explanation  of  these  fields. 


The  formats  for  assembly  station  costs  are  exactly  the  same  as  those 
t of  the  manufacturing  station  costs,  except  for  the  title  card.  For  each 

assembly  station,  one  needs  to  Input  the  cost  coordinates  only  once  although 
there  are  many  components  running  Into  the  station.  This  is  done  by  select- 
ing any  component  number  belonging  to  that  station  as  having  the  proper 
cost  function. 


Title  Card 


Col  U::ll 

\ 7 3 « 5 f.  ? t 

•••xr  \ii  ;•  3 ,.i'. /.■■"n  r rr  7."/-n  r \r,  .’i/  nr.- 

usmir  costs 

CuW'mI  •. 

Subsequent  Cards 

See  the  corresponding  format  of  MANUFACTURING  C0STS,  Section  4.3.5. 
4.4  Cost  Function  Input  Formats 

A cost  function  which  describes  the  cost  versus  the  quality  level 
is  required  for  every  raw  material,  manufacturing  station  and  assembly 
station.  Ten  types  of  codes  (0,1,..., 9)  are  provided  (See  also  Section 
6.2)  to  Identify  the  types  of  cost  functions.  Table  4.2  shows  the  types 


of  cost  functions. 

» 

Table 

4.2  Types  of  Cost  Functions 

Type 

Cost  Function 

l 

0 

Every  pair  of  cost  coordinates  is  individually 
given  to  the  program. 

1 

Standard  cost  function  supplied  by 

the  program 

2,. ..,9 

User's  cost  functions  in  functional 

form. 

U 

All  of  these  cost  functions  must  be  a monotonlcally  nonincreasing 
function  of  the  quality. 

In  Sections  4.3.4,  4.3.5  and  4.3.6  we  describe  the  input  formats  for 
raw  material,  manufacturing  station  and  assembly  station  costs.  All  of 
these  have  similar  formats  l.e.,  one  card  which  Identifies  the  component 
(and  stage)  and  the  type  of  cost  function;  and  subsequent  card(s)  where 
Interpretations  of  data  are  conditional  upon  the  type  of  cost  function 


* 

1 | 


specified  in  the  first  card.  The  interpretations  are  as  follows: 


For  a station  or  raw  material,  if  the  cost  function  is  of  type  0, 
then  the  number  of  coordinate  pairs  must  be  inputted  in  the  first  card. 
Subsequent  card(s)  will  contain  the  coordinates.  Five  pairs  of  coordi- 
nates are  allowed  on  each  card.  (See  Section  4.4.1  for  more.)  If  the 
cost  function  is  one  of  the  types  1,2,..., 9,  then  only  the  next  card  will 
be  considered  as  related  to  the  station.  The  cost  function  parameters  on 
this  card  will  be  used  by  SUBROUTINE  C0STFN  to  generate  cost  coordinates. 
(See  Section  4.4.2  and  4.4.3. ) 

4.4.1  Type  0 

When  the  cost  data  cannot  be  expressed  in  functional  form  or  only  a 
discrete  set  of  alternatives  are  to  be  considered,  the  known  coordinates 
are  Inputted  directly.  For  example,  if  for  Component  #3  of  the  raw  mater- 
ials, we  have  a choice  of  6 brands  of  the  same  kind  of  raw  material,  we 
can  represent  them  by  the  cost  function  of  Figure  4.5. 


Figure  4.5  Sample  Cost  of  Type  0 


25 


To  Input  data  In  the  RAW  MATERIAL  C0STS  section  for  this  component, 
the  first  card  Is: 


I i 

The  general  format  for  subsequent  card(s)  is: 

\ i 

r l 

$ represents  the  cost/unit  at  a specified  value  of  quality  level,  p. 

Hence  for  this  component,  2 cards  are  required  and  the  coordinate 
pairs  are  packed  left  justified  as  follows: 
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4.4.2  Type  1 

This  Is  the  standard  cost  function  supplied  by  the  program.  Once 
type  1 Is  specified,  the  next  card  has  the  following  format: 


tolian 

1 ? J 4 5 6 7 8 

1 2 345678901  ? 34  5678901  ? 34  5678901  ? 3456709G 1 334  567890  IT  34  56 ,890ir  34  56  78001^34  SCO0" 

10X 

rs.3 

fi.  3 

F5.J 

F5.3 

f$.3 

• 

b 

c 

*0) 

A(2) 

Comwnti 

where:  a * grid  size  of  p 

b * lower  bound  on  value  of  p(b>0) 
c * upper  bound  on  value  of  p(c>b) 

A(l),  A( 2)  * constants 

The  corresponding  equation  to  generate  the  cost  coordinates  Is: 
[$/un1ts]  * + A(2) 

Figure  4^6  Illustrates  the  cost  function. 
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4.4.3  Type  2, 3,. ..9 


If  the  user  wishes  to  use  a functional  form  other  than  Type  1 and 
does  not  want  to  put  all  the  cost  coordinates  in  as  for  type  0,  then  the 
user  may  modify  the  subroutine  C0STFN  to  specify  the  user's  functional 
form.  Section  6J?  describes  this.  The  user  may  specify  up  to  8 additional 
functional  forms  and  call  them  Types  2 through  9. 

If  one  of  these  types  is  specified,  only  one  card  containing  the 
values  of  parameters  for  that  functional  form  is  input.  Seven  parameters, 
A(l),  A(2),...,A(7),  are  allowed  to  be  used  in  any  of  these  cost  functions. 
The  format  is: 


C0l  U I 


12  3 4 5 6 7 8 

l?34r'C.78SQl?34507n0)2345r7890ir34r.67S?0123<567C?M234rX78‘;0123<!>67b<i01?:i:'..-?C'’'i 


10* 

■ 75. 3 

75. 3 

75.3 

F5.3 

75.3 

w 

75.3 

F5.3 

75.3 

75.3 

ft 

b 

c 

*0) 

A(2) 

A(3) 

A<4) 

A(5) 

A (6 ) 

A ( 7 ) 

Connects 

where  a = grid  size  of  p 
b = lower  bound  of  p 
c * upper  bound  of  p 

A(1 A(7)  * constants  to  be  used  in  SUBROUTINE  C0STFN. 

SUBROUTINE  C0STFN,  as  supplied  by  the  package,  has  the  functional 
equation: 

[$/unit]  = 1 

for  types  2,3 9 but  these  are  given  only  as  examples.  For  example, 

to  specify  a different  functional  form  for  Type  2,  the  statement  with 
statement  number  2 In  Subroutine  C0STFN  would  be  changed  to:  2 FN  * 

"your  functional  form  using  up  to  7 parameters  and  X is  the  quality  level". 
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5.0  OUTPUT 


This  chapter  provides  the  characteristics  of  the  outputs  as  well 
as  their  interpretations.  Two  options  of  the  outputs  are  given:  namely, 
the  complete  output  (intermediate  output  and  sumnary  output),  and  only 
the  summary  report  table.  Both  are  obtained  through  the  variable  NPRT 
in  the  main  program.  NPRT  equals  1 if  the  former  is  desired,  zero  if  the 
latter  is  desired.  The  default  option  is  complete  output.  Section  5.1 
gives  the  contents  of  the  summary  report  table  and  Section  5.2  gives  the 
contents  of  the  complete  output.  The  interpretation  of  the  complete  out- 
put and  the  diagram  of  the  optimal  solutions  to  an  example  problem  are 
discussed  in  Sections  5.3  and  5.4. 

5.1  Summary  Report  Table 

As  mentioned,  when  the  user  requests  only  the  final  solution,  the 
main  program  must  be  modified  to  make  the  variable  NPRT  = 0.  In  effect, 
the  summary  report  table  is  printed  immediately  after  all  input  images  and 
the  program  limitation  indicators  which  are  explained  in  Section  6.1.2. 

The  following  discussion  will  illustrate  the  output  for  the  example 
given  in  the  Appendix.  The  input  stream  and  the  complete  output  are 
presented  as  an  aid  to  the  user. 

The  summary  report  table  consists  of  sixteen  columns,  which  provide 
the  information  needed  for  each  station.  The  information  is  printed  sta- 
tion by  station  first  component  to  the  last,  starting  with  the  last  stage 
first  and  proceeding  to  the  first. 

The  following  gives  the  list  of  column  headings. 
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Table  5.1  Column  Headings  - Summary  Report 


1 234  56  7 89 

STATION  STAGE  C0MP0-  TYPE  SAMPLING  PLAN  PROCESS  ‘STATE  VARIABLE 


NENT 

OF 

STN 

N 

C 

FRAC-  -m- 

DEFECT 

A0QL 

: 

10 

11 

12 

13 

14 

15  16 

‘STATE 

VARIABLE 

STAGE 

MINIMUM 

SUB 

ALT  REMARKS 

| 

n 

— 

COST 

EXPECTED 

0PT 

SOLN 

i 

pi 

PT 

COST  F * 

IND 

IND 

Table  5.2  indicates  which  columns  are  used  for  each  type  of  station. 

Table  5.2  Output  Column  Usage  ii| 


Station  type  (x  - required) 

COLUMN 

NUMBER 

lllllil  nil——— 

, , .7/ 

First  line 
(PR0D  subr. ) 

the  following  lines 
(ASSMBY  subr.) 

1 

Station  name 

X 

X 

X 

2 

Stage  number 

X 

X 

X 

X 

3 

Component  number 

X 

X 

X 

X 

4 

Type  of  station 

A 

A 

i 

p 

5 

Sample  size 

X 

6 

Acceptance  number 

X 

7 

Process  fraction 
defective 

X 

X 

8 

Average  outgoing 
quality 

X 

9 

Average  outgoing 
quality  level 

X 

10 

Incoming  quality 

X 

X 

X 

X 

11 

Process  outgoing 
qual ity 

X 

X 

12 

Stage  cost 

X# 

X# 

X# 

13 

Minimum  expected 
Cost  (F‘) 

X 

X 

X 

14 

Suboptimal 

indicator 

X 

X 

X 

X 

15 

Alternative  solu- 
tion indicator 

X 

X 

16 

Remarks 

* 

* 
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**  There  are  two  types  of  information  for  each  assembly  station;  because 
for  every  "in- component"  assembly  station,  subroutine  ASSMBY,  (see  Chapter  3, 
package  overview  and  Appendix  1)  is  called  (m- 1 ) times  in  order  to  determine 
aggregate  qualities  then  subroutine  PR0D  is  finally  called.  Therefore, 
in  total,  m lines  of  information  are  recorded  for  each  assembly  station. 

# In  the  first  stage  and  for  any  type  of  station,  stage  cost  shown  ex- 
cludes the  raw  material  costs.  In  other  words,  the  raw  material  costs  can 
be  determined  by  subtracting  stage  cost  from  minimum  expected  cost  (F*). 

The  first  stage  of  a component  need  not  be  an  inspection  station. 

5.2  Complete  Output 

Besides  the  summary  report  table,  the  user  may  need  the  intermediate 
outputs  of  all  state  variable  values,  station  by  station,  for  all  stages 
as  well  as  three  tables  derived  from  the  INSPl  subroutine,  namely  the  ave- 
rage total  inspection  table,  sample  size  table,  and  acceptance  number  table. 

The  average  total  inspection  table,  the  sample  size  and  acceptance  number 
tables  are  entered  by  the  incoming  quality  and  the  average  outgoing  quality 
level  values.  The  information  will  serve  the  user  as  a source  for  the  con- 
ventional DP  iterative  solution  trace  back.  This,  in  turn,  enables  the 
user  to  obtain  the  solutions  manually  for  sensitivity  analysis. 

After  the  INSPl  output,  the  solution  table  of  each  station  is  printed, 
station  by  station,  component  by  component,  stage  by  scage,  all  starting 
from  the  first  to  the  last.  This  will  be  followed  by  the  summary  output 
as  described  in  Section  5.1 . An  inspection  station  or  manufacturing  station 
requires  one  table.  In  an  "m-component"  assembly  type  station,  two  types  of 
tables  are  printed  out.  The  first  ones  are  called  "aggregate  incoming  qua- 
lity" tables  and  (m-1)  tables  of  this  kind  are  printed.  The  last  table  contains 
the  same  kind  of  information  as  a manufacturing  station  table  and  refers  to 
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the  additional  defectives  generated  by  the  assembly.  Therefore,  for  each 
assembly  station,  m tables  are  printed.  All  stations  are  identified  by 
component  number,  stage  number,  station  type  and  station  name,  (if  any) 
in  headings.  Table  5_^3  gives  the  information  for  an  Inspection  station, 
as  conceptualized  in  Figure  2.1 . 


Table  5.3  Inspection  Station  Output 


COLUMN 

NUMBER 

VARIABLE 

DESCRIPTION 

REMARKS 

1 

AOQ 

Average  Outgoing  Quality 

AOQ  ^ AOQL 

State  variable 

2 

AOQL 

Average  Outgoing  Quality 

Limit 

3 

N 

Sample  size  "n" 

Decision  variable 

4 

C 

Acceptance  number  "c" 

Decision  variable 

5 

PIY 

Incoming  Quality  Level 
"Pj"  (Optimum  for  AOQ) 

State  variable  from 
previous  stage 

6 

STAGE  COST 

Inspection  and  Repair  Costs 

— 

Does  not  include  raw 
material  cost 

7 

F0 

Minimum  Average  Cost  Thus 

Far  (f„> 

f = f , + STAGE  COST 
n n- 1 

8 

ks  n 

Suboptimal  Indicator 

(see  Section  5.4) 

Table  5.4  gives  the  information  generated  for  a Manufacturing  station 
as  conceptualized  in  Figure  2.2. 

Table  5.5  shows  the  output  for  the  (m-1 ) tables  for  an  m-component 
assembly  station.  Figure  5.1  illustrates  the  notion  of  this  recursive 


aggregation  of  components. 


In  effect.  Figure  5.1  shows  the  computer  programming  aspect  which 

uar  rrncroH  in  Carf  i rvn  9 A Tha  flCCMRV  cuhr'Aiitl  no  i C ral  1 pH  f m—  1 \ timPQ  . 


according  to  Figure  2^3,  to  determine  recursively  the  minimum  cost  aggre- 
gate incoming  quality  for  that  assembly  station.  The  first  call  on  ASSMBY 
would  treat  Pj  1 as  PI1  and  Pj  2 as  PI2,  this  results  the  values  of  PAI 
with  the  values  of  optimum  PI1,  PIA1  and  optimum  PI2,  PIA2.  These  values 
are  stored  and  printed  in  the  first  table  (if  complete  output  is  desired). 
The  second  call  on  ASSMBY  would  treat  PAI  from  the  first  table  as  PI1  and 
Pj  ^ as  the  new  values  of  PAI  with  the  values  of  optimum  of  previous 
PAI,  PIA1  and  optimum  pj  PIA2.  The  calls  continue  in  this  manner  (m-1) 
times  where  the  PAI  from  the  (m-1)  st  call  represents  the  aggregate  incom- 
ing quality  for  that  assembly  station. 


Figure  5.1  Assembly  Station  Quality  Aggregation 


Table  5.4  Manufacturing  Station  Output 


COLUMN 

NUMBER 

VARIABLE 

DESCRIPTION 

REMARKS 

1 

PT 

Average  Outgoing  Quality  (py) 

State  variable 

2 

PIY 

Average  Incoming  Quality  (p  ) 
(Optimum  for  (py)  1 

State  variable  from 
previous  stage 

3 

PPY 

Average  Process  Quality  Level 
(P') 

Decision  variable 

4 

STAGE  COST 

Manufacturing  Cost 

Does  not  include  raw 
material  cost 

5 

F0 

Minimum  Average  Cost  Thus  Far 
<fn> 

f = f , + STAGE  COST 
n n- 1 

6 

LSUB 

Suboptimal  Indicator 

(see  Section  5.4) 
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Table  5.5  Assembly  Station  Quality  Aggregation  Output 


COLUMN 

NUMBER 

VARIABLE 

DESCRIPTION 

REMARKS 

1 

PA  I 

Aggregate  Incoming  Quality 
Level  (pAI) 

PAI  = PIA1  + PIA2 
- PIA1*PIA2 

2 

PIA1 

(aggregate)  incoming  quality 
(optimum  p^j  or  PI1) 

PIA1  = PI1  if  the 

1st  tab 
= PAI  of  pre- 
vious tab 
otherwise 

3 

PIA2 

Average  Incoming  Quality 
(optimum  PI2  for  PAI) 

State  variable  to  the 
previous  stage 

4 

FA  I 

Minimum  Aggregate  Cost 

FAI  = FI1  + FI2 

5 

ISUB 

Alternative  solution  indicator 

(see  Section  5.4) 

Table  5.6  shows  the  output  generated  by  optimization  with  respect 


to  p. 


Table  5.6  Assembly  Station  Decision  Variable  Output 


COLUMN 

NUMBER 

VARIABLE 

DESCRIPTION 

REMARKS 

i 

PT 

Average  Outgoing  Quality  (^) 

State  variable 

2 

PIY 

Average  Incoming  Quality  (pj) 
(optimum  for  PT) 

PIY  (PT)  = opt.  (PAI) 
PPY 

3 

PPY 

Average  Process  Quality  Level 
(p‘> 

Decision  variable 

4 

STAGE  COST 

Assembly  Cost 

Does  not  include  raw 
material  cost 

5 

F0 

Minimum  average  cost  thus  far 
<fn> 

f = fn  , + STAGE 
n "_1  COST 

6 

LSUB 

Suboptimal  Indicator 

(see  Section  5.4) 
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5.3  Interpretation  of  Outputs  and  An  Example 


This  section  will  discuss  the  interpretation  of  the  outputs  and  the 
manual  traceback  procedure.  The  example  of  Figure  5.2  will  be  used  as  an 
illustration.  Appendix  A contains  the  computer  output  for  this  problem. 

It  also  contains  the  input  card  images. 

When  the  solutions  for  a specified  final  AOQL  at  the  last  station 
(always  an  inspection  station)  is  required,  the  sunmary  report  table  would 
be  sufficient  for  the  source  of  actions  for  the  system.  The  cost  shown  on 
the  first  line  of  the  13th  column  represents  the  minimum  average  cost 
incurred  for  the  entire  system  for  that  particular  value  of  the  final  AOQL. 

Whenever  the  user  needs  to  know  the  average  total  inspection  or  desires 
to  examine  plans  for  final  AOQL  values  other  than  that  specified,  the  com- 
plete outputs  must  be  called.  The  latter  can  be  obtained  by  intermediate 
outputs,  station  by  station  and  the  basic  knowledge  of  dynamic  programming 
trace-back. 

5.3.1  Manual  Trace-Back  Procedure 

In  effect,  the  suimiary  report  table  can  be  obtained  manually  through 
the  intermediate  outputs,  station  by  station.  Fortunately,  for  a specified 
final  AOQL,  the  program  has  a certain  routine  to  perform  a trace-back 
internally.  In  what  follows,  we  will  demonstrate  how  the  user  can  obtain 
the  solutions  to  other  values  of  AOQL.  We  will  describe  the  procedure 
using  the  example  in  Appendix  A with  AOQL  = 0.050. 

Before  we  proceed  to  the  example,  the  procedure  is  recapitulated 
as  follow  : 

1.  It  is  advisable  to  have  the  structure  of  the  system  in  front  of 
the  user  to  reduce  ambiguities  in  connecting  the  stations  and  the  components. 
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2.  Start  from  the  last  stage  and  proceed  toward  the  first  stage. 

3.  Start  from  last  table  to  the  first  for  each  station. 

4.  Finish  all  stations  in  that  stage  before  proceeding. 

5.  Table  b_J_  below  summarizes  the  state  variables  used  to  connect 
stations. 


Table  5.7  State  Variables  Names  - By  Stations 


TYPE  OF  STATION 

OUTGOING 

STATE  VARIABLE 

INCOMING  STATE  VARIABLE 

REMARKS 

Inspection 

A0Q 

PIY 

Manufacturing 

PT 

PIY 

"n-component" 

PT 

PIY 

the  nth 
table 

Assembly 

PAI ( PI Y 

PAI(PIA1 

table) 

from  nth  table) 

from  (n-1  )st 

1 

1 

PIA1  and  PIA2* 

PIA1  and  PIA2* 

1 

1 

the  (n-1  )st 
table 

the  (n-2)nd 
table 

PAI (PIA1 

1 

from  2nd  table) 

l 

1 

PIA1 * and  PIA2* 

the  1st 
table 

* These  are  the  average  incoming  qualities  Pj  of  the  nth  component  to 

the  1st  component  entering  to  that  particualr  "n-component"  assembly  station, 

respectively. 

Figure  5.2  illustrates  the  structure  of  the  system  used,  for  an  example: 
The  final  AOQL  = 0.050. 

The  number  above  the  line  are  the  average  quality  values  (optimum  out- 
going and  incoming  qualities). 

Starting  from  the  last  table  of  the  intermediate  outputs,  component 
#4,  stage  #5  which  is  the  inspection  station.  The  specified  (system)  ave- 
rage outgoing  quality  limit  is  0.050  which  yields  the  A0Q  value  of  0.049, 
n = 71,  c = 6,  and  the  required  incoming  quality  (P I Y ) of  0.055. 
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COMPONENT 


STAGE  #1 


#2 


#3 


#5 


#4 


i .Qfia 


0.045 


2 .Q&5J 


3 -m 


4 .Q&Ll 


MfiSL 


.Q80J 


0-000 

: (Y\  ninn 

n nnn 

1 

S>-// 

Q>QQQ 


0.000 


\ 0.085 

o.ooo  L--^7 

f^\\0.055 

1 

n 

AOQ  = 

. 049-.A0QL 
^ .050 


Figure  5.2  System  Structure  - An  Example 


In  the  table  of  the  assembly  station  in  Stage  4,  the  state  variable 
used  is  PT  (the  average  process  outgoing  quality)  which  is  the  PIY  value 
obtained  previously,  i.e.  0.055.  Consequently,  the  required  average  in- 
coming quality  (PIY)  = 0.000  with  PPY  (the  average  process  quality  required) 
is  0.055. 

In  the  table  Component  5,  Stage  4,  the  state  variable  for  this  table 
is  PAI  (the  aggregate  incoming  quality)  = 0.000.  From  PAI  value,  the  PIA1 
is  0.000  and  the  PIA2  is  0.000  which  will  be  used  as  the  state  variable  to 
the  table  of  Component  5,  Stage  1. 

Entering  table  of  Component  4,  Stage  4 by  the  state  variable  PAI 
which  is  the  value  of  PIA1  = 0.000  obtained  previously,  the  PIA1  is  now 
0.000  and  PIA2  is  0.000  which  will  be  used  as  the  state  variable  to  the 
table  of  Component  4,  Stage  3. 
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Again  the  latest  PIA1  value,  i.e.  0.000,  is  used  for  the  state  var- 
iable PAI  of  the  table  of  Component  3,  Stage  4.  The  new  PIA1  is  0.000 
and  the  PIA2  is  0.000.  Since  the  present  table  is  the  first  one  for  this 
assembly  station  (Component  4,  Stage  4),  the  former  value  of  0.000  is  the 
state  variable  A0Q  value  for  the  inspection  station  of  Component  2,  Stage  3. 
The  later  value  of  0.000  is  the  state  variable  A0Q  value  for  the  inspection 
station  of  Component  3,  Stage  1. 

In  the  table  of  Component  4,  Stage  3,  an  inspection  station,  the  value 
of  state  variable  (AOQ)  is  0.000  obtained  from  the  PIA2  value  of  the  table 
of  Component  4,  Stage  4.  The  n and  c values  are  1000  and  zero  with  PIY 
value  (the  required  incoming  quality  level)  equals  0.085. 

In  the  table  of  Component  2,  Stage  3,  an  inspection  station,  the 
value  of  AOQ  is  0.000  obtained  from  the  PIA1  value  of  the  table  of  Com- 
ponent 4,  Stage  4.  The  plan  is  n = 1000,  c = 0,  and  the  incoming  quality 
level  is  0.100  (PIY). 

In  the  table  of  Component  4,  Stage  2,  a manufacturing  station,  the 
value  of  state  variable  Py  is  0.085  obtained  from  PIY  value  of  table  of 
Component  4,  Stage  3.  The  required  incoming  quality  level  (PIY)  = 0.000 
and  the  average  process  average  is  0.085. 

In  the  assembly  station  in  Stage  2,  entering  the  table  with  the  PIY 
value  obtained  from  the  table  of  Component  2,  Stage  3 to  the  state  vari- 
able Py  which  is  0.100.  The  required  incoming  quality  level  (PIY)  is 
0.045  and  the  process  quality  level  is  0.060. 

Using  the  value  of  PIY  obtained  just  above  (G  045)  as  the  state  vari- 
able PAI  of  the  table  of  Component  2,  Stage  2,  PIA1  and  PIA2  values  are 
0.045  and  0.000.  These  values  will  be  used  as  the  state  variable  values 
(AOQ)  for  the  tables  of  Component  1,  Stage  1 and  of  Component  2 Stage  1, 


respectively. 


* 


> 
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At  Stage  1,  all  stations  are  inspection  stations,  we  have  already 
obtained  the  values  of  state  variable  AOQ  for  tables  in  Component  5 up  to 
1,  which  are  0.000,  0.000.  0.000.  0.000,  and  0.045,  respectively.  The 
values  of  average  incoming  quality  for  raw  materials  are  0.065,  0.065, 
0.080,  0.065  and  0.080,  respectively. 

Figure  5^_3  gives  the  solution  values  to  the  system  structure. 


Raw 
Mat'  1 


Figure  5.3  Solution  to  the  Example 
5.4  Suboptimal  Conditions  and  Alternative  Solutions 

At  any  station  there  exists  a suboptimal  indicator  variable  for  each 
value  of  the  resulting  quality.  The  indicator  is  0,  if,  at  that  average 
outgoing  quality,  the  value  of  decision  variable  is  a global  optimum;  1 
if  it  is  a "possible"  local  optimum.  The  variable  is  KSUB  for  Inspection 
and  LSUB  for  manufacturing  or  assembly.  For  an  inspection  station  the  sub- 
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optimal  condition  is  caused  when  the  average  incoming  quality,  Pj,  called 
for  by  the  algorithm  (see  Chapter  2)  is  greater  than  the  maximum  average 
incoming  quality,  Pj  max,  for  that  particular  stationl  This  means  that 
the  global  optimum  P*  might  exist  beyond  the  PT  value,  had  the  maximum 
average  incoming  quality  been  greater  than  the  current  P.  . 

Another  type  of  suboptimal  indicator  (ISUB),  again  zero  or  one,  is 
involved  in  assembly  stations  where  the  ASSMBY  subroutine  is  called  for. 

In  effect,  the  subroutine  computes  the  aggregate  incoming  quality  from  a 
pair  of  components  coming  into  assembly  stations.  ISUB  equals  1 if  an 
aggregate  incoming  quality  formed  by  two  sets  of  average  incoming  quality 
is  greater  than  the  maximum  values  of  average  incoming  quality  in  either 
set,  i.e. 

PAI  > min  [max  P 1 1 , max  P 12] 

Furthermore,  the  subroutine  gives  the  indications  of  alternative  solu- 
tions in  the  summary  report  table.  The  number  indicates  how  many  alterna- 
tive combinations  of  Pj^  and  P^  values  the  particular  PAI  has  that  yield 
the  same  incoming  costs. 


I 


1 (see  Ch.  3) 


I max 


t 


I raw 
AOQL  x 


FMAX 


: No  previous  station 
: otherwise 
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6.0  POSSIBLE  MODIFICATIONS  TO  THE  PROGRAM 


The  program  package  as  supplied  can  be  easily  modified  to  accomodate 
bigger  size  problems  and/or  other  cost  functions.  Larger  problems  are 
discussed  in  Section  6.1.  Other  functional  forms  for  cost  functions  are 
discussed  in  Section  4.4.3.  If  the  input  data  is  on  cards,  see  Section 
6.2. 

6.1  How  to  Accomodate  a Larger  Problem 

By  changing  the  dimensions  of  arrays  and  vectors  in  the  MAIN  program 
one  can  increase  the  size  of  the  problem  that  the  program  can  solve.  Twelve 
parameters  are  used  to  control  the  amount  of  storage  area  required  in  the 

computer.  They  are  indicated  and  defined  in  Table  6.1. 

The  listing  of  the  MAIN  program  is  in  Appendix  B.  By  altering  the 
DIMENSI0N  Statements  in  the  MAIN  program  according  to  the  appropriate 
parameters,  the  size  of  the  problem  that  the  program  can  handle  is  altered. 

It  is  advised  that  for  a large  problem,  the  program  should  be  run, 

first,  by  adjusting  only  the  MAXC0M  and  MAXSTG  parameters  to  fit  the  pro- 
big*.  Mhe  program  will  determine  the  proper  sizes  of  the  other  parameters 

for  adjustment  in  the  next  run.  The  reason  behind  this  approach  is  because 
the  task  of  estimating  the  sizes  of  the  various  parameters  is  a rather 
time-consuming  process  and  the  allocation  of  too  big  a storage  may  be  too 
costly. 
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Table  6.1  Parameter  Definitions 


Parameter 

MAXC0M 

MAXSTG 

01 

J2 

J3 


J5 

J6 

J7 

J8 

J9 

J10 


Definition 

Maximum  number  of  components 

Maximum  number  of  stages 

Total  # of  cost  coordinate  pairs  for  all 
RAW  MATERIAL  C0STS 

Maximum  number  of  cost  coordinate  pairs 
for  a station  or  a raw  material  component 

Total  # of  cost  coordinate  pairs  for  pro- 
cessing stations  and  assembly  stations 

Total  # of  inspection  stations 

Number  of  values  of  A0QL 
(FINAL  A0QL)  * (FMAX)  , 

State  Variable  Grid  Size 

Number  of  unique  Pj  values  used  to  generate 

tables  of  ATI  N and  C.  Note:  J6  < < J7 

J1  + J5 

(#  of  components  in  the  problem)  *J5 

/Maximum  ft  of  components  goingv*J5 
\into  an  assembly  station  / 

Maximum  of  (J2,  J6) 


Other  parameters  that  may  need  to  be  adjusted  for  input/output  devices 
are  in  Table  6.2. 


ft 


Table  6.2  Input/Output  Device  Parameters 

Parameter  Definition 

NRD  Input  1/0-unit  capable  of  rewinding 

NWR  Output  1/0-unit 

NDS  I/O  unit  for  random  access  file 

NPRT  fl  full  output  desired 

10  partial  output  desired 


6.2  Rewinding  of  Input  Data  File 


As  shown  in  Figure  4J_  (Schematic  Diagram  of  Input/Output  files). 
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it  is  required  that  the  input  data  for  a problem  be  put  on  a sequential 
data  file  which  can  be  rewound  by  the  program.  The  purpose  of  this  fea- 
ture is  for  the  program  to  copy  the  data  file  image  on  the  output  for 
user  verification,  (see  the  example  in  Appendix  A).  The  program  rewinds 
the  file  and  reads  it  again  to  interpret  the  data.  It  is  possible  to  elimi- 
nate the  rewinding  step  if  one  will  sacrifice  the  printing  of  the  input 
data  image. 

If  the  input  data  is  on  cards,  the  job  control  language  associated 
with  this  program  will  be  simplified  if  the  rewinding  step  is  eliminated. 

In  order  to  accomplish  this,  statements  with  sequence  numbers  IPT  00039 
to  IPT  00050  (inclusive)  of  Subroutine  INPUT  must  be  deleted. 


Note  1 The  program  as  distributed  will  solve  a system  with  up  to  14 
components  and  15  stages. 

Note  2 It  took  less  than  one  minute  of  cpu  plus  1/0  time  to  solve  the 
example  problem  in  Appendix  A. 
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